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In the present paper we examine gravitational wave asteroseismology relations for /- 
modes of rapidly rotating neutron stars. An approach different to the previous studies is 
employed - first, the moment of inertia is used instead of the stellar radius, and second, the 
normalization of the oscillation frequencies and damping times is different. It was shown 
that in the non-rotating case this can lead to a much stronger equation of state independence 
and our goal is to generalize the static relations to the rapidly rotating case and values of 
the spherical mode number l > 2. We employ realistic equations of state that cover a 
very large range of stiffness in order to check better the universality of the relations. At 
the end we explore the inverse problem, i.e. obtain the neutron star parameters from the 
observed gravitational frequencies and damping times. It turns out that with this new 
set of relations we can solve the inverse problem with a very good accuracy using three 
frequencies that was not possible in the previous studies where one needs also the damping 
times. The asteroseismology relations are also particularly good for the massive rapidly 
rotating models that are subject to secular instabilities. 

PACS numbers: 


I. INTRODUCTION 

A direct detection of gravitational waves is a major goal in physics for several decades. Nat¬ 
urally the investigation of the strong sources of gravitational radiation is being done in parallel. 
Different dynamical processes connected to neutron stars, such as neutron star mergers and ro¬ 
tational instabilities, are considered as promising sources. The final goal is to use the observed 
gravitational wave signal in order to extract the characteristic parameters of the observed objects, 
such as mass, radius, rotational rate, etc., and eventually to constrain the nuclear matter equation 
of state (EoS) 00 ]- 

In the present paper we will concentrate on the neutron star oscillations with an emphasis on 
the region of the parameter space where the Chandrasekhar-Friedman-Schutz (CFS) instability 
can develop, i.e. certain nonaxisymmetric modes can become unstable due to the emission of 
gravitational waves. This is one of the most promising scenario for the emission of strong gravita¬ 
tional radiation where the signal can reach even above the Advance FIGO sensitivity |[a0]. More 
specifically we will reexamine the gravitational wave asteroseismology with rapidly rotating neu¬ 
tron stars by considering a different parametrization of the EoS independent relations. 

Studying gravitational wave asteroseismology of oscillating neutron stars originates with the 
papers [3—5] where the mode oscillation frequencies and damping times of static neutron stars 
were related to their mass and radius. It was found that the asteroseismology relations do not 
depend on the EoS up to a large extend and the inverse problem can be solved with a very good 
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accuracy. Later these results were extended by including other modern realistic EoS 0. The 
generalization of the relations to the case of rapidly rotating neutron stars and larger values of 
the spherical mode number l was done recently in 0 , 0 ]- These results extend up to the Kepler 
limit which means that asteroseismology can be performed for a much larger class of objects and 
more specifically for neutron stars that are subject to the CFS instability. In a similar fashion the 
frequencies of the emitted gravitational wave signal during binary neutron star mergers can be 
connected to the properties of the post-merger supramassive neutron stars 00]. 

IQ ] connect the oscillation mode frequency and damping 


All the relations obtained in 


time to the neutron star mass, radius and rotational rate. But a different approach was proposed 
in Ell which steamed out of an earlier study by Lattimer and Schutz Ell - to use the moment 
of inertia instead of the stellar radius. In this way the obtained relations are much more EoS 
insensitive compared to the previous results. Similar relations were also obtained in Ell . But the 
results in [11, [13j] are limited to the nonrotating case. Our main goal in this paper is to extend 
them to the rapidly rotating case in order to see if the EoS universality is preserved and whether 
this approach is more beneficial compared to the standard case 0,0]. 


II. ASTEROSEISMOLOGY USING THE MOMENT OF INERTIA 


Our aim is to present the results as a brief report. For this reason we will not go into details 
about the method of calculation of the mode frequencies and damping times, but instead we will 
give only the most important points. We refer the reader to jz] where a detail presentation of the 
methodology can be found (see also 00,0). 

We obtain the oscillation modes by performing a time evolution of the linearized perturbation 
equations in the Cowling approximation, where the perturbations of the metric are neglected. The 
reason for these simplifications is that the full nonlinear simulations are extremely computer and 
manpower demanding 0 and solving the linearized perturbation equations without any addi¬ 
tional approximation is not done until now. But we should note that even though this approxima¬ 
tion can lead to relatively large deviations of the frequencies (in the range 10 — 30% depending on 
the compactness and the mode number /) it was shown in Istlzt] that it gives not only qualitatively 
but in some cases also quantitatively good results for the asteroseismology relations. 

The numerical code we use for the time evolution of the perturbation equations was developed 
in 0[0. The background neutron star solutions are obtained with the RNS code 117 ] that has 
been proven to be reliable for rapid rotation. The mode frequencies are easily obtained after 
performing a Fourier transform on the computed time series. Obtaining the damping time on 
the other hand is more involved. Since we are working in the Cowling approximation, where the 
background metric is not perturbed, the damping (or the growth in the case of CFS instability) of 
the modes due to the emission of gravitational waves can not be calculated directly 1 . Instead one 
can apply approximate Newtonian formula, where the emission of gravitational waves is related 
to the mutlipole moments of the neutron star. It was shown that at least in the nonrotating case 


1 It is questionable whether this damping/growth time can be accurately calculated using a time evolution code in 
general, since the inevitable numerical dissipations in many cases are much stronger than the dissipation of energy 
due to gravitational wave emission. 
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this formula gives good results for the damping time compared to the exact general relativistic 
values llil . 

We use three EoS that cover a very large range of masses and radii - EoS A 0, WFF2 0] 
and L 101 . The first one is very soft reaching a maximum mass of only 1.67 Mq. The second 
one is a standard modern realistic EoS that fulfills all the observational constraints on the neutron 
star mass and radius iEtI E^I . EoS L on the other hand is a very stiff (one of the stiffest available 
in the literature) with M max = 2.72 Mq and typical radii around 14 — 15 km. The set of EoS is 
chosen to cover such large range of stiffness for a reason - we want to verify up to what extend 
our results are insensitive to the EoS. For practical purposes one can of course use a narrower set 
of EoS that falls into the preferred range of masses and radii according to the observations. This 
can significantly reduce the error in the fitting formulae. Since we are calculating the modes in the 
Cowling approximation, we can not reduce the errors of our asteroseismology formulae below the 
deviations coming from this approximation. That is why we decided not to give separate fittings 
for a set of EoSs that fall into the preferred range of neutron star masses and radii. Our purpose 
here is more to offer an alternative way for doing asteroseismology of rapidly rotating neutron 
stars and to demonstrate its validity over a large range of EoS stiffness. 

Throughout the paper we will use the following dimensions for the different quantities. The 
mass M is measured in solar masses Mg, the frequencies (the angular rotational frequency O as 
well as the mode oscillation frequency a) are measured in kHz, the damping/growth time - in 
second and the moment of inertia I is defined as I — I cgs / 10 45 g cm 2 , where I C g S is the moment of 
inertia in cgs units. In this way for example the quantity r], that is widely used in the paper and 
roughly speaking has a meaning of an effective compactness of the star, takes the form 


V = 


M 


(. 


i. 


cgs 


Mq J \ 10 45 g cm 2 


(1) 


A. Oscillation frequencies 

Here we will present the results for the asteroseismology relations of the oscillation frequen¬ 
cies. In the previous studies of asteroseismology with rapidly rotating neutron stars 0,0 two 
step relations were constructed that map the oscillation frequencies of rotating neutron stars to 
their mass, radius and rotational rate. As a first step we have relations between the normalized 
oscillation frequencies and the rotational rate derived for different values of Z, where the oscilla¬ 
tion frequencies are normalized to the corresponding values in the nonrotating limit. As a second 
step the nonrotating frequencies are expressed as a function of the neutron star average density 
similar to fl]. 

Here we follow a different approach. We are going to use the moment of inertia instead of 
the stellar radius and a different normalization similar to [11]. It was shown that this makes 
the asteroseismology relations much more insensitive to the EoS. First, we will focus on counter 
rotating /-modes with l = m that are potentially CFS unstable. In Fig. [I] we plot the Z = m = 2 
/-mode oscillation frequencies measured in inertial frame of reference 2 a‘ ,nsl normalized to the 


2 The inertial frame frequency a\ is connected to the comoving frame frequency cr c by the relation cr c = v; + mCl/ln, 
where m is the azimuthal mode number. 
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FIG. 1: The potentially unstable l — m — 2 /-mode oscillation frequencies measured in inertial frame as 
functions of the parameter //. 


neutron star mass M, as a function of the parameter //. Different colors correspond to sequences 
with fixed values of the normalized rotational parameter Q = Mfi and different styles of the 
symbols correspond to different EOS. The data range from the nonrotating limit Q = 0 to the 
rapidly rotating case with Q = 25 (in units [M© kHz]). In the latter case the models have rotational 
frequencies in the range / = 1.2 — 1.5kHz and they do not have a stable nonrotating limit. 

The first thing one can notice is that even though the chosen EoS are very different, the data 
points for a fixed Q fall very well on a single line. Moreover the slopes of the fitting lines vary by 
a small amount for different values of O. It turns out that the dependences can be well approxi¬ 
mated by a fitting formula of the type 

a“ nst M = (ci] + a 2 tl + a^Ci 2 ) + ( b\ + b 2 tl + b^Cl 2 )!/. (2) 

The coefficients cii,..,bj, are given in Table 1 for different values of Z. 

TABLE I: The coefficients a\,..,b 3 in the fitting formula ||2), for the potentially unstable branches with l — 
in — 2,3 and 4. 



«i 

a 2 

a3 

bi 

b 2 

b3 

/ 

= m 

= 2 

-1.76 

-0.143 

-6.65 x 10~ 3 

3.64 

-4.36 x 10~ 2 

2.00 x 10’ 3 

z 

= m 

= 3 

-2.55 

-0.236 

-1.63 x 10” 2 

4.67 

-6.63 x 10~ 2 

5.79 x 10” 3 

z 

= m 

= 4 

-3.21 

-0.365 

-2.47 x 10~ 2 

5.55 

-7.33 x 10~ 2 

9.09 x 10~ 3 


An interesting fact, we should note, is that the coefficients a\,..,b 3 can be expressed as linear 
functions of the mode number / with a good accuracy (only the error in the fitting of b 2 is slightly 
larger). Thus we obtain 

af nst M = [(-0.332 - 0.725/) + (0.085 - 0.111 l)Cl + (0.0112 - 0.00903/)G 2 ] + 

[(1.755 + 0.955 /) + (-0.0165 - 0.0149 Z)Q + (-0.00501 + 0.00355 Z)Q 2 ] ij. (3) 
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An important consequence of this dependence is that with a good accuracy one can say that the 
relations given by eq. (2|) for different l are linearly dependent. This reflects the possible ways of 
solving the inverse problem that will be discussed in the next section. 

The above given asteroseismology relations have several advantages compared to the ones 
considered in [6,@]. First, they are much simpler and the rotational frequency Q enters at second 
order compared to the third order relations given in 10, 7], Also the deviations from EoS univer¬ 
sality are smaller 3 . But there are of course also disadvantages compared to the original approach 
IQ, I7j]. Here the asteroseismology relations are prone to the deviations coming from the Cowl¬ 
ing approximation, whereas in the two step relations in it was shown that the normalized 
relations between the oscillation frequencies and the rotational rate are very close to the results 
coming from the full nonlinear GR calculations 00. 

We will present here also briefly results for the asteroseismology relations for the co-rotating 
stable modes with l = —m. Part of our motivation comes from the fact that by employing a specific 
normalization of the quantities, we can solve the inverse problem by using three frequencies (both 
stable and unstable). This was not possible using the relations in 03 where in order to obtain 
the mass, radius and rotational rate independently one has to use at least one damping time of a 
mode. 

The relations for the stable branches are similar to the potentially unstable case. But it turns 
out that for solving the inverse problem it is better to use the oscillation frequencies in comoving 
frame u c instead of o Also the dependence on O can be well approximated with a linear function 
compared to the second order polynomial in the potentially unstable case. The normalized de¬ 
pendences between the l = — m = 2 oscillation frequencies in comoving frame and the parameter 
ij are shown in Fig. [2] for different values of Q. 

We use a fit very similar to the one for the counter-rotating modes given by eq. |2]l with the 
only difference that the dependence on Q is linear: 


cj^}I sA — (tq T / 72 O) T (£q + f?2fd)/y. 


(4) 


The coefficients a.\, 02 , b\ and £>2 are given in Table II. We present results only for the l = —m = 
2 modes, and not for higher values of /, since in general the modes with lower value of l are 
better emitters of gravitational waves. For the counter-rotating potentially unstable modes we 
considered also the / > 2 case because the CFS instability might develop easier for higher l 0,0, 
14]. But we do not have such a mechanism for the co-rotating modes and it is generally accepted 
that only the co-rotating stable / = 2 modes will emit significant amount of gravitational waves. 


TABLE II: The coefficients a\, a? ar| d b in eq. Q for the stable branch with l = — m = 2. 



a\ ci2 bi &2 

1 — —m — 2 

-1.66 -0.249 3.66 0.0633 


3 As a matter of fact the normalized relations between the oscillation frequencies and the rotational rate in are 
also quite independent of the EoS. But the second set of relations, namely between the nonrotating frequencies and 
the neutron star average density, can lead to larger deviations from the fits as discussed in 0,113]. 
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FIG. 2: The stable l — —m — 2 /-mode oscillation frequencies in comoving frame as functions of the 
parameter )]. 


B. Damping times 


Deriving asteroseismology relations for the damping time t^w is not as straightforward as for 
the oscillation frequencies for the following reasons. The damping time can be expressed as 


tgw = 


1 dE 
2E~dt' 


(5) 


where E is the mode energy and the energy carried by gravitational waves dE / dt is proportional 
to the squared l + 1 time derivative of the corresponding multipole moment. Applying the New¬ 
tonian multipole formula for the gravitational wave emission one can show that lb. 18, 2dh 


1/tgw °c cf Z+ V C/ 


( 6 ) 


where <7/ and <J C are the frequencies in the inertial and corotating frames respectively. As one can 
notice tgw —> 00 when <7, —> 0. This means that tgw can not be expressed directly as a function of 
the stellar mass, moment of inertia and rotational rate similar to eq. © because then the peculiar 
behavior of Tgw around the <7, = 0 point (i.e. around the transition from stable to CFS unstable 
regime) can not be reflected. Instead Tgw can be expressed as a function of (7/ similar to the 
previous studies |7,14]. In order to obtain relations insensitive to the EoS one has to of course use 
a proper normalization of the quantities, that is explained and justified below. 

Making some approximate back of the envelope calculations on the basis of eq. (0, and taking 
into account that the parameter ;/ can play the role of an effective compactness instead of ( M/R ), 
one can obtain 


V 



cx Mcr. 


(7) 


Even though this formula is based on very simplified assumptions, it offers a good way for nor¬ 
malizing the damping time asteroseismology relations. 
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FIG. 3: The normalized damping time ?/ (M/t// 2 ) ^ ] as a function of the normalized oscillation frequency 
Mcr“ nst for the l = m — 2 and l — m — 4 /-modes. 


The dependences of the normalized /-mode damping time as a function of the normalized 
oscillation frequency for the potentially CFS unstable branches with l = m = 2 and l = m = 4 are 
shown in Fig. [3j One can easily notice that the dependences pass through the origin of the axes 
where both a\ mst = 0 and 1 /tgw = 0 according to eq. ©. Therefore we can approximate the data 
with a second order polynomial of the form 


r] = Cl {Ma? nst ) + c 2 (M<7"" sf ) 2 , 


( 8 ) 


where C\ and c 2 are constants. If we substitute ( Ma \ inst ) from eq. Q into eq. we will finally 
have the desired relation between the damping time from one side and the neutron star mass, mo¬ 
ment of inertia and rotational rate from the other. As one can see from Fig. [3j the approximation 
is very good for l = m = 4 but the errors in the l = m = 2 case can be larger. But these large errors 
are observed mainly in the stable part of the graph (with tgw > 0)- That is why we have chosen 
to give separately the values of Ci and c 2 obtained from the fit of the unstable part only. Moreover 
this is the regime where the CFS instability is operating and can offer a promising mechanism for 
the emission of detectable gravitational waves signal 0,0. The values of the coefficients C\ and 
c 2 for / = in = 2, 3 and 4 are given in Table III both for the full sequence of data and for the CFS 
unstable part of the branches where Tqw < 0 and u, < 0. 

Let us now compare the relations derived here with the previous studies. A single dependence 
between the normalized damping time and the normalized rotational frequency was derived in 
0,0 for the potentially unstable branches with l = m = 2,3 and 4. Thus it was in a way more 
general than the corresponding relations in the present paper given by eq. ©. But as we ex¬ 
plained above, the relations in 00 are two step, and the normalized dependences have to be 
supplemented with relations between the static damping times and the stellar mass and radius. 
These static relations are different for different values of l and the deviations from EoS universal¬ 
ity can reach large values. That is why if we sum up the advantages and disadvantages in both 
case, we can conclude that the relations given here are at least as good as the ones in 0 , 0 . Of 
course in order to determined more precisely which set of relations is better one has to examine 









8 


TABLE III: The coefficients c\ and c 2 in eq. © for different values of l. 

Full sequence: 



Cl c 2 

1 = m — 2 

l — m — 3 

l — m — 4 

0.644 2.07 x 10~ 2 

0.246 2.68 x lO^ 3 

0.147 1.91 x 10~ 4 


Unstable part: 



Cl c 2 

l = m — 2 

l — m — 3 

l — m — 4 

0.403 -7.64 x 10 -2 

0.217 -4.82 x 10” 3 

0.144 -3.15 x 10” 4 


also the inverse problem and see how accurately we can derive the stellar properties from the 
observed mode frequencies and damping times. This will be done in the following section. 

We should note that the Cowling approximation can introduce very large deviations in the 
damping times. For example if the frequency increases with say 30%, the damping time decreases 
by a factor of three (for the l = 2 mode) due to the strong dependence of t gw on a, according to 
eq. ©. That is why relaxing the Cowling approximation is very important for the exact calcula¬ 
tion of the damping times and the related gravitational wave asteroseismology. Such a project is 
underway. 


III. SOLVING THE INVERSE PROBLEM 

In this section we will consider the inverse problem, i.e. obtaining the stellar parameters using 
the observed oscillation frequencies and damping times. The previous results for rapid rotation 
showed that using the normalizations in 0,0 one can not obtain the mass, radius and rotational 
rate of the star independently using only the observed oscillations frequencies due to the specific 
form of the relations used there. Instead one has to use also the damping times. But observing the 
damping times is much more difficult compared to the oscillation frequencies and the errors are 
expected to be considerably bigger. 

The normalization presented in this paper allows us to improve the picture a little bit. Un¬ 
fortunately one still can not solve the inverse problem using three frequencies of CFS unstable 
modes with different /. The reason is that as we said above a\, .., b 3 in eq. CD can be expressed 
with a very good accuracy as linear functions of / which makes the corresponding coefficients for 
/ = 2,3 and 4 linearly dependent. Instead one can supplement the unstable 1 = m mode frequen¬ 
cies with a frequency of a stable corotating mode having m < 0. 

In Table IV we present the results after solving the inverse problem in four different potential 
scenarios: 

1. We are able to observe the l = m = 2,3 potentially unstable mode frequencies together with 
the l = —m = 2 stable mode frequency. 

2. We are able to observe the l = m = 2,4 potentially unstable mode frequencies together with 
the l = — m = 2 stable mode frequency. 
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3. We are able to observe the l = m = 2,3 potentially unstable mode frequencies together with 
the / = m = 2 damping time. 

4. We are able to observe the l = m = 2,4 potentially unstable mode frequencies together with 
the l = m = 2 damping time. 


We will represent the accuracy of the inversed asteroseismology relation on two neutron star 
models that are massive and rapidly rotating and are thus subject to the CFS instability. Such 
models can be the outcome of core-collapse or binary neutron star merger and they are promising 
sources of gravitational waves BB- Moreover it was shown that exactly the rapidly rotating 
massive models are the biggest challenge to the asteroseismology relations presented in BB and 
lead to large errors when solving the inverse problem. In the cases 3. and 4. where damping time 
is used as input information, we employ the relations for the unstable part of the asteroseismology 
relation given in Table III. 


TABLE IV: Solutions of the inverse problem using either three frequencies or two frequencies and a damp¬ 
ing time. Two different neutron star models with EoS WFF2 that are subject to the CFS instability are 
chosen. The rounded percent deviations from the exact values are given in brackets. 



M[M 0 ] 

I/10 45 [g/cm 3 ] 

n[kHz] 

M[M 0 ] 

l/10 45 [g/cm 3 ] 

n[kHz] 

Exact 

2.02 

2.80 

8.66 

2.34 

3.04 

9.62 

& £ 7f!g & afL 2 

2.02 (0.2) 

2.82 (1) 

8.63 (0.3) 

2.39 (2) 

3.17(4) 

9.48 (1) 

of "f & c& ofl 2 

2.02 (0.2) 

2.82 (1) 

8.63 (0.3) 

2.37(1) 

3.17(4) 

9.49 (1) 

of =2 & of™* & 

2.15 (7) 

2.41 (14) 

9.70 (12) 

2.35 (1) 

2.72 (10) 

10.12 (5) 

t rfl 2 & of™* & rf™* 

2.17(7) 

2.41 (14) 

9.75 (13) 

2.34 (0) 

2.72 (10) 

10.09 (5) 


From the data in Table IV one can make the following conclusions. First, the neutron star 
parameters can be obtained with a good accuracy using the given combinations of potential ob¬ 
servables. This is particularly true for the first two cases with three observed frequencies where 
the errors are very small. Therefore the relations presented here are more accurate than the ones 
in BB especially in the massive rapidly rotating neutron star case. Also the advantage we have 
here is that the stellar parameters can be obtained using three frequencies. This is important since 
the damping times will be much more difficult to be obtained from the detected signal and the 
observational errors are supposed to be much bigger compared to the oscillation frequencies. An¬ 
other important point is that the relations presented here are good also for very massive rapidly 
rotating models that were the biggest challenge to the relations in B B- Such neutron stars are 
of particular interest since they are subject to the CFS instability and are supposed to emit strong 
gravitational radiation in the early stages of their life BB- That is why obtaining asteroseis¬ 
mology relations that can be used for solving the inverse problem with a good accuracy for such 
models is important. 

It is natural to expect that the accuracy when solving the inverse problem can be improved in 
two ways. First, as we commented above, is to use a restricted set of EoS, because in this paper 
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we employed a quite broad set in order to test the EoS independence of our relations and some 
of the EoS are even not in agreement with the current observational constraints. Another way of 
improving the results is to make a second step in solving the inverse problem explained in detail in 
laid]- Namely, we can use the obtained values of M, Q and ?/ and rederive all the asteroseismology 
relations for values of the parameters close to these particular M, Q and //. Afterwards we can use 
the new relations and obtained updated values of M, Q and rj. This procedure can be repeated 
more than once. 

To conclude the section let us discuss a possible modification of the asteroseismology relations 
presented here. A key parameter in the relations is rj which plays the role of generalized compact¬ 
ness. We used it since it was shown that it leads to the largest known degree of EOS independence 
Q. But one can obtain similar asteroseismology relations if the true compactness M/R is used 
instead 4 . The way of building the relations will be practically the same as the one used in the 
paper, only ;/ has to be substituted with M/R. Our rough calculations show that the obtained 
relations in the case of M/R are also EOS independent up to a large extend. We will not give the 
explicit relations and numbers here, since the main idea of the paper is not to give as complete set 
of relations as possible, but instead to demonstrate an alternative way of doing asteroseismology 
and test its accuracy. We plan to derive a more complete set of relations in the future when we 
drop the Cowling approximation. 


IV. CONCLUSIONS 


In the present paper we considered a new way of parameterizing the gravitational wave as¬ 
teroseismology relations in the case of rapidly rotating neutron stars. We concentrated on the 
/-modes that are one of the most efficient emitters of gravitational radiation. Our study is moti¬ 
vated by the results in where it was shown that if one normalizes the oscillation frequencies 
and damping times in a certain way and uses the moment of inertia of the star, the asteroseismol¬ 
ogy relations become much more insensitive to the nuclear matter equation of state compared to 
the commonly used relations involving the neutrons star average density |4]. The results in 111] 
are limited to the static case and here we extended them to rapid rotation and also to values of 
the spherical mode number l > 2. This is particularly important since the /-modes of neutron 
stars can become CFS unstable for rotational frequencies above roughly 80% of the Kepler limit 
QfiSJ thus becoming efficient emitters of gravitational radiation. The case of 1 > 2 on the other 
hand is important since the standard mass neutron stars (with mass say M < 2.OM0) have much 
larger CFS instability window for modes with l = m = 3 and l = m = 4 compared to / = m = 2. 

In our studies we used three different equations of state that cover a very large range a stiffness 
- from very soft EOS with small maximum mass, to very stiff EoS with large maximum mass and 
large radii. Even though some of these EOS are not in the preferred range of neutron star masses 
and radii according to the observations, they are very useful for testing the EOS independence 
of our results. It turns our that after choosing an appropriate normalization, the obtained aster¬ 
oseismology relations are indeed EoS independent up to a large extend for both the oscillation 


4 The asteroseismology relations using the compactness M/R were shown to be quite EOS independent in the static 
limit 
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frequencies and damping times. 

There are several important differences with the previous results on rapidly rotating neutron 
stars obtained in 0,0 . The main one is that the asteroseismology relations in [6,0 are two step. 
First, we have relations between the normalized oscillations frequency (or damping time) and the 
normalized rotational frequency, where the normalizations involve the oscillations frequency and 
damping times in the nonrotating limit. Second, we have relations connecting the oscillations 
frequency and damping times in the nonrotating limit to the neutron star mass and radius. The 
normalization considered in the current paper is different - we have single relation connecting the 
oscillation frequency (or damping time) to the neutron star parameters (mass, moment of inertia 
and rotational frequency) for every value of 1. The relations in the present paper also look simpler 
- they are either linear of quadratic with respect to Q compared to the second and third order 
relations in 0 , 0 . Another difference is that we are using the moment of inertia instead of the 
stellar radius, since it was shown in [11] that this leads to better EOS independence. On the other 
hand one can well constrain the nuclear matter EOS using the moment of inertia 1121. 

The same methodology as the one described in the present paper can be used to derive aster¬ 
oseismology relations that employ the stellar radius, as discussed in the previous section. Our 
rough calculations show that these relations are also quite EOS independent (one just has to sub¬ 
stitute the effective compactness // with M/R). 

At the end we concentrated on different ways to solve the inverse problem, i.e. obtain the 
stellar parameters from the observed oscillation frequencies and damping times. The most im¬ 
portant thing to point out is that using the relations presented in the current paper, we were 
able to solve the inverse problem with a very good accuracy using three oscillations frequency, 
that was not possible in the previous studies 0,0 where the use of both the mode frequencies 
and damping/growth times was required. This is quite beneficial, since extracting the damp¬ 
ing/ growth times from the gravitational wave signal might be very difficult and the accuracy 
will be much lower compared to the oscillation frequencies. When solving the inverse problem 
we concentrated mainly on massive models that are rotating rapidly, i.e. the primary candidates 
for the CFS instability. The derived relations are particularly good in such cases that represent 
the biggest challenge to the relations in 0,0 , where the errors increase significantly for massive 
rapidly rotating models. 
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